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ON THE CONTINUATION OF A POTENTIAL GAS 
FLOW ACROSS THE SONIC LINE 


By Lipman Bers 


SUMMARY 


Simple, sufficient conditions are given under which a two-dimensional 
steady compressible flow can be continued across the line of Mach number 1 
as a continuous supersonic flow. Methods for the actual computation of 
the flow are described. The problem is of importance in the theory of 
transonic flows past airfoils. 


INTRODUCTION 


The theory of mixed subsonic and supersonic flows of a compressible 
■fluid is as yet very incomplete. As a contribution to such a theory, 
the problem of continuing a two-dimensional subsonic gas flow across a 
line along which the speed of the flow reaches that of sound is consid- 
ered in this paper. Rather simple, sufficient conditions under which 
the continuation is possible and is uniquely determined are given 
( section l) . The problem is of interest in connection with transonic 
flows past obstacles (see section 7)» 

The mathematical nucleus of the problem is the solution of a Cauchy 
problem for a partial differential equation of mixed elliptico-hyperbolic 
type. This problem is solved in the appendix by a method which is 
believed to be intrinsically simpler than the previous treatment by 
Christianovitch in a paper unavailable in this country^ and by Frankl 
(reference l) . The theoretical discussion contained in the appendix 
(sections 5 and 6) yields methods for the actual computation of the flow 
in the supersonic region. 

This investigation was carried put at Syracuse University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


■''Christianovitch, S. A.: On Supersonic Gas Flows. CAHI Rep. 

No. 5^3, 1941. (in Russian.) 
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SYMBOLS 


a 

Fa 

K 


K* 

% 

M 

n 

q 

^cr 

^max 

s 

S 

t 


^max 


The following notation is used in the ANALYSIS: 
speed of sound 

solution of hypergeometric equation (27) 

coefficient in Chaplygin's equation, defined hy 
equation (20) 

step function approximating function K 
function defined hy equation ( 38 ) 

Mach number 

normal to the curve T pointing toward the domain 
speed 

critical speed 
m axim um value of q 

arc length measured on F 
length of T 

function of speed (see equation (l8)) 
value of t for M = » 


A 
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u,v 


x,y 

X(s),Y(s) 

7 

r 

A 

e 

p 

p cr 

T 

T cr 

f 

V 

cu 


components of velocity vector 
coordinates in physical plane 
coordinates of points on T 
exponent in pressure— density relation 
curve in. physical plane along which M = 1 
image of r in holograph plane 
domain in which a subsonic flow is given 
angle between velocity vector and x— axis 
density 

critical density 

variable defined by equation (28) 
value of t for M = 1 
velocity potential 
stream function 

angle between tangent to r and x— axis 
used over symbol, value on r 


0 


used as subscript, value for q = 0 
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The following notation is used in the appendix: 


end points of segment along which Cauchy data are given 


AjBjTjNjT' ,N' ,T",N" constants defined in section 3 of the appendix 


D(a/b) 


domain defined in section 1 of the appendix 


coefficient in equation (Al) 


functions approximating function K 


complex— valued functions 


independent variables 


y (n)^y(n) functions defined by equations (Al8) 


complex variables 


i -z(n) "formal powers" defined in section 4 of the appendix 


constants defined by equations (A3^-) 


matrix defined by equation (A9) 


initial values of ■\jr and its no rmal , derivative 


function connected with by equation (AT) 


dependent variable in equation (Al) 
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ANALYSIS 

1. The Continuation Theorem 


The main result of this investigation reads as follows: 

Given a domain A in the x, y-plane, whose "boundary contains 
a continuously curved arc T (see fig. l); given a subsonic 
potential gas flow in A; and if 

(a) the Mach number of the flow on F is equal to unity, 

(b) the components of the velocity vector and their first—' 
order derivatives are continuous on F, 

(c) the normal derivative of the speed on F is different 
from zero, and 

(d) no streamline of the flow cuts F at a right angle; 
then the flow can be continued across any subarc of r as a 
potential supersonic flow without weak discontinuities and this 
continuation is uniquely determined. 

This theorem contains no general statement concerning the extent 
of the supersonic region obtained; however, in each special case the 
method of continuation described as follows determines the extent. 

By a compressible potential flow without weak discontinuities is 
meant a flow described by two twice continuously differentiable 
functions 0(x, y) and ty(x, y) (potential and stream function) 
satisfying the differential equations 


M = fo 

Sx p By 

M. M 

Sy P Sr 


( 1 ) 
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where Pq is the stagnation density and p the density given hy 


P = PQ 





1 


(2) 


Here 7 is the exponent in the pressure— density relation; 8 q is the 
speed of sound at a stagnation point; and q, the speed, is 


q. 



(3) 


It is well known that in the supersonic region flows with weak 
discontinuities (discontinuities in the second derivatives of the 
velocity potential) are possible. In fact, there is no good physical 
reason to assume that in a supersonic flow the velocity- components 


S0 

= ®o r" 

dx 

► 


v 




possess partial derivatives. (The continuity equation and the condition 
for the existence of a potential may be expressed in terms of integral 
relations.) Under precisely which conditions flows with weak disconti- 
nuities can be obtained by continuing subsonic flows across the sonic 
line seems to be a very important but rather delicate question; it is 
intimately tied up with the question of the necessity of conditions (b) 
to ( d) . 

If condition (c) is violated along an arc of T, then the velocity 
vector is constant along this arc. This follows from equation (8) . It 
seems very probable that in such a case the velocity vector must be 
constant everywhere. Whether the continuation theorem remains true if 
condition (c) is violated at isolated points remains an open question. 
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Condition (d) is probably essential for the uniqueness of the 
continuation. If at some point P of r the streamline Is normal 
to r, then the Mach lines at this point are tangential to r , for 
at F the Mach angle is 90°. Such a situation arises, in a symmetrical 
Laval nozzle if the transition from subsonic to supersonic speeds occurs 
along a line extending across the nozzle (flow of Mayer type, see 
fig. 2). In a remarkable paper (reference 2) Fr ank! investigated this 
case and showed that the continuation of the subsonic flow into the 
supersonic region is not uniquely determined. He showed also that in 
this case the supersonic flow may exhibit weak discontinuities. 


2. Reduction to 
In order to prove 


a Cauchy Problem in the Holograph Plane 
the theorem, let the equation of the curve 


r be 


x = X(s) 
y = Y(s) ” 
0 < s < 


(5) 


where s is the arc length of P measured in the direction of 
traversing T with the domain A to the left. By hypothesis X"(s) 
and Y"(s) exist and are continuous, and X'2 + y'2 _ it follows 
that the angle 


0 )(s) = tan 1 y {? . ) ■ 
X«(s) 


between the tangent to r at x = X(s), y=Y(s) and the positive 
x-direction is a continuously differentiable function of s. (See 
fig. 1.) 

Set 


■10 


u — iv = qe' 


(7) 
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The value of any function F(x,y) at x = X(s) , y = Y(s) will he 
denoted by F(s). Thus 0(s) = 0[x(s), Y(s)J. Nikolskii and Taganov 
(reference 3) showed that 


0 *( 8 ) 



8n 


(8) 


where 8/8n denotes differentiation in the direction of the normal 
to r pointing toward the domain A. (This formula follows easily 
from the equation of motion written with • q and 0 as- independent 
variables.) Since M = 1 on r and M < 1 in A (M=q/a being 
the Mach number) it follows that 8q/8n < 0, so that by hypotheses (c) 
and ( d) 


0*(s) < 0 (9) 


It is no loss of generality to assume that 


2* > 0(0) > 0(S) > 0 


( 10 ) 


If this condition were not satisfied, T could be divided into a 
finite number of overlapping arcs along each of which 0 would change 
by less than 2rt. Inequalities ( 9 ) and (10) imply that F has a one- 
to-one image in the holograph plane, which is the circular arc F^ 
given by 


<L = % t , 9(S)^0<0(O) 


( 11 ) 


where q cr , the critical speed, is 


tr 
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By the definition of the potential and stream function. 




= r- (u to + v dy) 

a o 


(-T to + n dy)J 


( 12 ) 


On V, 


(V ry 

u = q cr cos ® 

<v 

v = q cr sin 6 
to = ds cos cd 


dy = ds sin cd 



so that 



These equations show that ^(s) ^and i{r(s) are twice continuously 
differentiable. By hypothesis 9(a) is a continuously differentiable 
function. By virtue of inequality (9) the inverse function 
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s = s*(e) 


exists and hence is continuously differentiable. It follows that 


0*(e) = ?[s*(0)]l 


► 

^*( 0 ) =t[s*(0)] 

A» 

are continuously differentiable functions of 0,. 0(s ) < 0 < 0(0) • 

Assume now that the flow in A has been continued across F . 
Since 


(14) 


d(u,v) = ^5(^,0) 
5(x,y) 5(x,y) 


it follows that on T 


Mill = Jhf cob 2 (9 - «,) ( 15) 

5(x,y) \5n/ 

In fact, on r 

5(q,0) _ 80 _ 50 

5(x,y) 5s 5n 5n 5s 

and 5q/5s = 0, whereas 50 /5s is given by equation (8). By 
hypotheses (c) and (d) expression (15) cannot vanish. Hence 


< o 

d(x,y) 


for all points sufficiently close to F. It follows that a sufficiently 
small neighborhood of any subarc of T has a one— to— one image in the 
hodograph plane, that is, a neighborhood U of a subarc of In the 
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holograph plane the functions f and ijr satisfy the veil— hnovn 
Chaplygin equations 


80 . IPO 
80 P 84 


M 1 - M 2 PO 8t 

8g q p 80 ^ 


( 16 ) 


vhile on they satisfy the initial conditions 


0(4cr,e) = 0*(9) 

► 

*(4cr,e) = **(e) y 


(IT) 


for 0(S) < 0 < 0(0). 


3 . Existence and Uniqueness of the Solution 
in the Holograph Plane 

In the preceding section the continuation of a subsonic flov 
across the sonic line has been reduced to the initial— value problem, 
equations (l6) and (17) . Introduce the new independent variable 



( 18 ) 


(Note that q = 0 corresponds to t = — q = q cr to t = 0, 
and q = q max = a^ \j2/(y — l) to some finite positive value of t.) 
In the 0,t— plane, 0 and satisfy the equations 
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M M 

► 

M = _K(t)^i 

St 8s j 

■where 

E(t) . k=*L 
2 


(19) 


(20) 


The initial— value problem, equations (l6) and (17) , goes over into 
a Cauchy problem for equations (19) • .The theory of this Cauchy problem 
is given in the appendix. From the theorems stated in section 3- of the 
appendix it follows that the initial— value problem, equations (l 6) 
and ( 17 ) , possesses a unique solution in the domain bounded by the 
arc and two characteristics of equations ( 16 ) . As is well known, 

these characteristics are epicycloids. In this domain (see fig. 3) , 
which is called the "characteristic triangle determined by 

the functions 0 and i|r will be continuously differentiable, provided 
that 0*(0) and T]r*(0) have this property. 

Remark .— Mote that if ijr (considered as a function of 9 and t) 
is twice continuously differentiable, then equations (19) imply that 
tJt satisfies Chaplygin’s equation 


K(t)£-£ + =• 0 

& 2 St 2 


'( 21 ) 


In this case the initial— value problem considered may be formulated as 
a Cauchy problem for equation (21) : ' 


*(9,0) = t(0) 

/ 


%(e,o) = v(e) 
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where 

t( 9 ) = f*(e) 
v(e) = 0*’(e) 


In general, however, need not he a twice continuously 
differentiable function in- the supersonic part of the hodograph plane. 


4. Existence and Uniqueness of the Solution 
in the Physical Plane 

In order to complete the proof of the continuation theorem, it 
must he shown that the functions j ' and obtained in the preceding 
section can he transferred (in a unique way) to the physical plane 
(more precisely, to a part of the physical plane adjacent to r and 
lying outside A) and that in the x, y-plane these functions possess 
continuous derivatives of the second order which assume the "correct" 
values along f . 

The functions 0 and obtained from the Cauchy problem coincide 
along P with the potential and stream function of the flow given 
in A and transferred to the hodograph plane by means of the mapping 


u = u(x,y) 

► 


v = v(x,y) 


( 22 ) 


It follows from equations (12) that the inverse mapping is given by 


x + 


l7 ’ !i °J s rfe+ "t 4 *) 


(23) 
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or. 


s + iy = ao / — 



+ 1 h) 


d0 



Since the functions 0 and i|r have been continued into the 
supersonic part of the hodograph plane , this integral is meaningful 
even there. Moreover, it is path— independent. Under the assumption 
that the partial derivatives }qq j $Qq> • • • exist and are continuous, 

this follows immediately from Chaplygin's equations (l 6). But, in the 
case considered, this assumption cannot he made. Consider, however, the 
mapping from the hodograph plane to the 0,^— plane. By virtue of 
equations (l6) 


8 ( 0 ,*) 
8(0, q.) 


= -£~ 
Pq<L 


/Mf + 

G£ 

/Mf 

\30/ 

p 0 (l - M 2 ) 

\8q/ 


It follows (for instance, from relations (9) and (13)) that this 
Jacobian is different from zero in some neighborhood of the arc T h . 

Hence the mapping from the hodograph to the 0,*— plane is locally one- 
to-one, so that it is sufficient to establish the path-independence of 
the integral, equation (23), in. the 0, ''Jr — plane . This leads at once 
to the conditions 



d_ 

dq 



= 1 - M 2 


Using the relation 
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these conditions may be transformed into 


ci© _ P Bq 
W ~ P 0 4 

► 

5© PoC 1 ~ ^0 ciq 

cty pq 50 


( 24 ) 


Interchanging the dependent and independent variables, it is seen that 
system (2k) is equivalent to equations (l 6). Thus equations (16) insure 
the pathr-independence of the integral, equation (23), even without any 
assumption about the existence of second-order derivatives of 0 
and ijr. 


Thus equation (23) defines a mapping also in the supersonic part 
of the holograph plane. The Jacobian of this mapping is easily 
computed to be equal to 



so that for q = c^. 


MZj7) 

8(u,v) 



(25) 


But by equations (14) 


mue) = 32 isl - (. g V 72 

ds d© \y + 1 j 


cos(9 - m) , 
S*(s) f 


so that 


hiisxl d 0 

8(u,v) 
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on the axe defined by equation (ll) and hence also in some neighborhood 
of any subarc of this arc. It follows, that such a neighborhood is 
mapped onto a neighborhood of a subarc of F in a one— to— one manner. 
The functions 0 and obtained by solving the initial— value problem, 
equations (16) and (17), can now be transferred to the x,y-plane by 
means of the mapping defined by equation (23). 

The derivatives of 0 and in the x, y— plane can be computed 

by using equations (l6) and (23). (See reference 4.) This direct 
though lengthy computation yields exactly equations (12). Thus 0 
and satisfy equations (l). Finally, equation (23) shows that x 
and y are continuously differentiable functions of u and v, so 
that u and v are continuously differentiable functions of x 
and y. Thus 0 and possess continuous partial derivatives of the 
second order in the x, y-plahe. The proof is now complete. 


5. Computation of the Solution 

The theorems proved in the appendix implicitly contain methods for 
the effective conputation of the solution. Thus theorem 2 shows that 
the solution of the Cauchy problem may be represented as in infinite 
series of particular solutions of the Chaplygin type. (The idea of the 
following proof is already contained in a remark by Frankl concerning 
the Tricomi problem; see reference 5-) 

Chaplygin (reference 6) showed that a particular solution of the 
second— order equation obtained from equations (l6) by eliminating 0 
can be found by setting (for any positive constant a) 


♦ - (26) 


where q max is the maximum speed possible for the given gas: 



and F a (T) is any solution of the hypergeometric equation 
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t (1 - T)F a "(r) + [l + a - (a a + b a + l) T ] F a *(T) - a^Fjr) = 0 (27) 


the constants and h a being determined by the equations 


a„ + b„ - a — 

a a y x 


a-h- = - S&L+il 

2( r -l) 


It, is seen that the potential corresponding to equation (2 6 ) is 
given by 


0 ^^^CPo/^C^/lmax) 11 §’a( T ) ~ (2T/a)F a ' (t)J 


where 



(28) 


Let FojCt) and F a jj(r) denote the particular integrals of 
equation (27) which for 


t = T 


Icr 


cr 


-max 


satisfy the initial conditions 
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F 


all 


= 1 


W 


(^J 


Set 


B aj =T“/% aJ (T) 




M,n J 


then 


fej - 

- ieiaS V T > , 


(29) 


( 30 ) 


axe pairs of solutions of equations (l6) satisfying the initial 
conditions 
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^ Vi = 0 for 4 3 = 1 


aj 


( 31 ) 


0 a j = °* Vi = 1 for ^ = 5 cxj J = 11 


■aj 


Consider now the Cauchy problem, equations (l6) and (17), ass umin g 
for the sake of simplicity that 


9(S) = -T 
9(0) = T 

0< T < it 


Set 



(32) 


and expand the initial data 0*(9) and ty*(9) in Fourier series 


0*(e) 



n=0 


^( 0 ) =y~ c n cos 0^9 + d n sin a n 0 f 
n=0 


( 33 ) 


T < 9 < T j 


From the way the particular solutions, equations (30), have "been 
constructed it follows that the series 
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CO 

0 = ao + ^> (&n cos 0^0 + "bn. sin o^A^itr) + 

n=l 


> (c n sin 0^0 - d n cos a^A^utr) 
n=l 


> 


♦ 



,(c n cos On© + d n sin a n0) B a n Il( T ) + 


(- 0 ^ sin oCjj0 + "bn cos c^Bo^iCt) 

n=l 


(34) 


represent a formal solution of the initial-value problem considered. 

Theorem 2 of the appendix not only shows that the series converge 
absolutely and uniformly (within the "characteristic triangle") to the 
desired solution but also permits an estimate of the error committed by 
replacing the infini te series by finite partial sums. Theorems 4 and 5 
of the appendix contain similar statements concerning the convergence 
of the differentiated series. 

Remark 1.- While the theorems in the appendix refer to an 
equation of the form of equation (2l) , it is clear that these theorems 
imply similar ones for system (l6) . It may be noted that the convergence 
of series (34) was proved by a general method; no special properties of 
the hypergeometric functions were used. 

Remark 2 .— It should be noted that if the arc r h on which the 
initial values are given exceeds 



then the epicycloids drawn from the extremities of the arc do not 
intersect. In this case the characteristic triangle must be replaced 
by the domain bounded by the initial arc, the two epicycloids, and an 
arc of the circle q = q max . (See fig. 4.) 
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6. Computation of the Solution by a Second Method 

The expansion in a Chaplygin series may not be the most efficient 
way of computing the solution of the initial— value problem in the 
hodograph plane, since it requires the preliminary computation of the 
functions Aq, and B a , and because of the possibility of slow 
convergence. Another method is suggested by theorems 3 and 6 of the 
appendix. 

Let the equation satisfied by the stream function be transformed 
into the form of equation (23). Divide the interval 0 < t < t^y 
into subintervals 


t y t <C ty 

v = 1, 2 . . . n 
0 = to < ti < . . . < t n - t mp ^- 

and set 


(35) 


K*(t) = min K(t‘) 


(36) 


ty_l — t < ty , ty_]_ < t * < ty 

The equation 


K*(t)^ 00 + tfct = 0 


(37) 


is an approximation to equation (21). Since K*(t) is a piecewise 
constant function ( step function ) , it is equivalent to the one- 
dimensional wave equation in each strip ty < t < t v+ j. A solution of 
equation (37) will be required to possess continuous partial derivatives 
of the first order. The potential 0 corresponding to a given 
solution of equation (37) will be continuous but its first derivatives 
will possess discontinuities. The characteristics of equation (37) are' 
given by the equations 
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-K*(t) 


dt = Constant 


9 + 




dt = Constant 


which represent polygonal arcs in the 0,t— plane. (See figs. 5 and 6.) 

For equation ( 37 ) the solution of the Cauchy problem, 
equations ( 2.K ) , is a trivial matter. This solution maybe constructed 
by repeated use of the classical D'Alembert formula (as described in 
detail in section 6 of the appendix) . Another possibility is to apply 
the method of the preceding section. 

Theorems 3 and 6 assert that for a sufficiently fine subdivision, 
equations ( 35) , the solution of equation (37) will be an arbitrarily 
close approximation to the solution of equations (2^) with the same 
initial values. 

Instead of using the approximating step function, equation ( 36 ) , 
one may use the function 


K*(t) = max K(t*) ( 38 ) 

t v _l ^ ^ t v> V-l ^ t * £ *v 

In this case K*(t) s 0 for 0 < t < tj, so that in this strip i[r is 
a linear function of t. 

In some cases it might be preferable to use as the approximating 
coefficient a piecewise linear function. In this case the 
approximating equation will be of the type of the Darboux— Tricomi 
equation 


^xx + V = ° (39) 

in each strip- ty < t < t v+ y. The solution of the Cauchy problem for 

the Darboux— Tricomi equation can be expressed either by integral 
formulas (see reference 7 ) or by infinite series involving Bessel 
functions . 
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7. Application to Transonic Flows past Obstacles 

The theorem stated in section 1 is of a purely local character. 

While in the holograph plane the flow can be continued throughout the 
whole characteristic triangle, the transition to the physical plane may 
be possible only in the immediate neighborhood of the sonic line r. 

This may be due either to the existence of limiting lines or simply to 
the fact that a simply covered domain in the holograph plane corre- 
sponds to a multiply covered domain in the physical plane. 

In some important cases, however, it can be asserted that the 
continuation of the subsonic flow by the method described in sections 2 
to 5 yields the whole supersonic region which is of relevance to the 
problem considered. 

Consider a transonic flow past a closed body with a subsonic stream 
Mach number (and without shock waves). The supersonic regions Sp 
and S 2 (see fig. 7) are bounded by solid walls and transition lines 
along which M = 1. For a supersonic region thus bounded, Nikolskii 
and Taganov (reference 3 ) proved an important theorem stating that the 
region S possesses a one— to— one image in the holograph plane. 

In this proof, Nikolskii and Taganov make use of the remark that 
all Mach lines in S must possess points in common with the transition 
line ( the line M = l) . While they give no formal proof of this fact, 
a proof can be supplied without difficulty. The domain S is simply 
covered by a family of smooth curves, the streamlines of the flow, and 
the solid wall bounding S belongs to this family. A Mach line having 
no points in common with the transition line would have to originate 
and terminate at the solid wall. A simple argument, essentially equiva- 
lent to Nolle ' s theorem, shows that such a Mach line would be tangent 
to a streamline at at least one point P. But this is impossible 
since at P the Mach angle would have to be 0, which corresponds 
to M = 00, p = 0. 

Consider now the hodograph image of S. It is bounded by the 

arc Fh of the circle q = q cr (corresponding to the transition line 
in the physical plane) and the curve W^ (the hodograph image of the 
solid boundary) . Assume that is not contained within the character- 
istic triangle determined by which is bounded by and two 

characteristics of the hodograph equations ( i .e . , two Busemann 
epicycloids) . Then at least one of these Busemann epicycloids 
intersects Wft at two points, and contains an arc of an 

epicycloid which has no points in common with F^. This, however, is 
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impossible since tbe epicycloids are tbe holograph images of the Mach 
lines. Thus is contained within the characteristic triangle. By 

virtue of the previous results , this implies the following uniqueness 
theorem: 


Two transonic (subsonic at infinity) flows past an obstacle 

are identical if they are identical in the subsonic region. 

Similar theorems are true for flows past curved walls and for 
Taylor type flows in nozzles (figs. 8 and 9)* Iu the case of a Mayer 
type nozzle flow, however, the subsonic flow does not determine the 
supersonic flow uniquely. 

The preceding remarks are of interest in connection with various 
attempts to construct transonic flows by the "correspondence method." 

The correspondence method was first used by Chaplygin for jet 
problems (where, however, it is identical with the solution of the so- 
called direct problem) . The correspondence method may be expressed in 
various analytical forms j however, except for the very original method 
of Bergman (reference 8) the basic procedure is always the same, as is 
clearly pointed out by Gelbart (reference 9) • It consists of associ- 
ating with a solution of the C auchy^-Ri emann equations (representing an 
incompressible flow) a solution of Chaplygin’s hodograph equations. 

In some cases the resulting compressible flow is of the same general 
character as the initial incompressible flow. To date, the most 
detailed application of this method to flows past obstacles is that 
given by Tsien and Kuo (reference 10) . The continuation of a sub- 
sonic flow across the sonic line invariably occurs in the application 
of the correspondence methods to transonic flows. 

While the mathematical problems connected with the correspondence 
method are as yet largely unsolved, the results of the present investi- 
gation show that once the solution of Chaplygin's equation in the 
subsonic part of the hodograph plane is determined, the continuation 
of the solution in the supersonic region is uniquely determined and can 
be effectively computed. If this continuation does no’t yield the 
desired result in the physical plane, the choice of the solution in the 
subsonic domain must be abandoned. 


CONCLUDING REMARKS 


The continuation of a given subsonic flow into the supersonic 
region has been discussed. The corresponding problem in the subsonic 
region is of little interest, since in the subsonic region the 
differential equations are of elliptic type with analytic coefficients 
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and have analytic solutions . Thus if the velocity potential 0 is 
known up to a line r, the continuation of 0 across r is a 
problem in analytic continuation. It seems improbable that any 
interesting and general statements can be made concerning this cane. 

In the supersonic region the equations are of hyperbolic type and the 
continuation of a given flow across a line r reduces to the stan dar d 
Cauchy problem. The theory of this problem is known (see reference 11 , 
p . 326 , and the literature quoted therein) . It would be very 
interesting, however/ to obtain results concerning a problem converse to 
the one considered here} that is, assume that a supersonic flow is 
known up to the line T on which the Mach number is equal to 1, and 
determine under which conditions the flow can be continued into the 
subsonic (elliptic) region. 


Syracuse University 

Syracuse, N. Y., September 1, 19^7 
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APPENDIX 

THE CAUCHY PROBLEM FOR CHAPLYGIN'S EQUATION 


This appendix contains "the proofs of the theorems used to 
establish the preceding results . For the convenience of the reader, the 
notation standard in the theory of partial differential equations is 
used. Note that x and y no longer refer to the physical plane. 
Equation (Al) is essentially the same as equation (2l) . 


1. Introduction 

The aim of this investigation is to integrate the differential 
equation 




+ ^ = ° 


(Al) 


where 

K(0) = (f 

K(y) < 0 ► 
for y < 0 

under the initial conditions 

tKx.,0) =_t(x) 
i}r y (x,0) = v(x) ► 
a 5L x < b 


(A2) 


' (A3) 
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t and V being given functions. The solution is required in the 
domain D(a,b) consisting of the segment a < x < b of the 
axis y = 0 and of the part of the half plane y < 0 bounded by two 
characteristics of equation (Al) passing through the points (a.O) 
and (b,0) . 

Chaplygin’s equation for the stream function as a function of the 
hodograph variables can be brought into the form of equation (Al), the 
domain y< 0 corresponding to the supersonic region. In the following 
considerations, however, no use will be made of the special form of the 
function K. 

The classical theory fails to treat the Cauchy problem, 
equations (Al) to (A3), except for analytic t and V, because the 
initial data are given along the line y = 0 where the equation ceases 
to be of hyperbolic type. For one special case, that of the Darboux- 
Tricomi equation 


^XX + *yy = ° (A4) 

the solution has been known for decades. In a recent paper, 
unavailable in this country, Christianovitch treated the case of the 
Chaplygin equation. According to Frankl (reference 1) he had to 
assume t(x) and v(x) to be analytic and had to use all their 
derivatives in order to obtain the solution. Frankl (reference 1) 
considered the case of nonanalytical Cauchy data for the equation 


ytacx + ^yy + a ( x >y) 1 i f x + + c(x,y)t = 0 (A5) 

assuming the coefficients a, b, and c to be analytic. This 
treatment is based on special properties of the Riemann function for 
the Darboux— Tricomi equation and of the hypergeometric function and 
does not seem to lend itself easily to numerical computations. Note 
that if K(y) is analytic and K'(0) ^ 0, equation (Al) can be 
brought to the form of equation (A5) by introducing 



, as a new independent variable. 
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The present treatment, while limited to equations of the form 
of equation ( Al) , is of an entirely elementary character, requires only 
very simple regularity conditions for the function K, and remains 
valid even where a reduction to the form of equation (A5) is impossible 

(e.g. , for K = yl/3) . Besides establishing the existence, uniqueness, 
continuity, and differentiability of the solution, it will be shown 
that the solution depends continuously on the Cauchy data (as noticed 
by Frankl) , as well as on the coefficient K. These continuity 
properties imply two effective methods for the approximate computation 
'of the solution: Either by representing it as a series of particular 

solutions of the Chaplygin type or by replacing the coefficient K by 
a piecewise constant function, in which case the determination of the 
solution becomes quite trivial. In fact, the approximation of K(y) 
by step functions constitutes the main tool of this investigation. 


2. Statement of the Problem 

The function K(y) will be assumed to satisfy the following 
conditions : 

(1) K(y) is defined and nondecreasing for y 5: 0 

(2) K(y) is continuous at y = 0 

(3) K( y) is negative, piecewise continuous, and possesses a 
piecewise continuous derivative for y < 0 

Since the only interesting case is K(o) = 0, the condition 
that K(y) be nondecreasing involves no serious loss of generality, 
(it is satisfied in the case of Chaplygin's equation.) 

A function t|t( x,y) will be called a regular solution of equa- 
tion (Al) if ijr possesses continuous partial derivatives of the first 
order throughout its domain of definition, and piecewise continuous 
partial derivatives of the -second order satisfying equation (Al) . 

The Cauchy problem mentioned in the Introduction may be stated as 
follows : 


Problem A . Given two twice continuously differentiable 
functions t(x) and v(x) for a < x < b. To determine a 
regular solution ijr(x, y) of equation (Al) defined in D(a,b) 
and satisfying the initial conditions , equations (A3) • 

It will be convenient to consider instead a more general problem 
(suggested by gas dynamical applications). A regular solution ilr of 
equation (Al) determines (except for an additive constant) a 
function cp( x, y) such that 
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<P 2 = *y 


cp y = -K(y)t 2 


(AT) 


(Note that the derivative cpy might possess discontinuities.) The 
complex— valued function 


v(z) = cp + iy 

> 

z A x + i y 


is called a Z-monogenic function, where 


1 1 

Z = Z(k) = 

1 K(y) 


(A8) 


(A9) 


Equations (A7) imply 


(A10) 


for every closed rectifiable curve T contractible to a point. It will 
be convenient to call every continuous complex— valued function, 
equations (A8) , satisfying equations (A10) a Z-monogenic function. 

Part of the theory of Z-monogenic functions (references 12 and 13) 
(originally developed primarily in the elliptic case K > 0) remains 


denotes the coefficient-matrix of equations (AT) 
that 


J cp dx — 


Ky dy = 0 


> 


y dx + cp dy = 0 
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valid with this new definition. In particular, every £— monogenic 
function w = <p + ity possesses a E— integral 

w(z) = f Z v(£). d s £ = r <P *= - & d y + 1 / ♦ + <P dy (All) 

J Zq « Z 0 V Z 0 

¥(z) is again a E-monogenic function, and 


w = 


5 ¥ 

Sx 


(A12) 


Now, if i]r is a solution of problem A, then (for a proper- 
normalization of qp) 


qp(x, 0 )= (v(!) d|, a <x <b 

Ja 


(A13) 


This remark suggests the following problem: 

Problem B . Given two continuous functions t(x) and v(x) 
for a < x < b. To determine a E-monogenic function w = cp + iijr 
defined in D(a,b) and satisfying the initial condition 

v(x) = J v(|) + 1t(x), a<x<b (Al 4 ) 

It is easy to see that if w is the solution of problem B, and 
if i}r possesses continuous first and piecewise continuous second 
derivatives, then ^ is a solution of problem A. 


3- Statement of the Results 

The main results of this investigation are: 

Theorem 1 (Existence and uniqueness of the solution) . 
t(x) and v(x) be arbitrarily given continuous functions 


Let 
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defined for a ^ x < b. There exists one and only one 
Z(K)-monogenic function v = Cp + ity defined in D(a,b) and 
satisfying the Cauchy condition 


. -w(x) 



V(t) d| + iT(x), 


a < x < h 


(A15) 


Theorem 2 (Continuous dependence on the initial values) . 

The function w = cp + it defined in theorem 1 satisfies in D(a,b) 
the inequalities 


|cp| < AT + BN 

► 

|t | £ T + |yj N 


(A16) 


where 


T = max |t(x)| 

N = max |v(x)J 
a < x < h 
A = A(y) = 

B = B(x,y) = x - a + |y|A(y) 


Theorem 3 (Continuous dependence on the coefficient of the 
equation) . Let K n ( y) he a sequence of functions satisfying 
conditions (l) to ( 3 ). Let w n = + if n he the S(K n )-Hnonogenic 

function satisfying condition (A15), t and v being fixed 
continuous functions. If 
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lim Ea(y) = K(y) 

n->oo 


uniformly for 0 > y > -H, where H is the least upper hound 
of |y| in D(a,b), then 


lim w n (z) = w(z) 
n — 


uniformly in D(a,b), w being the E(K)— monogenic function 
satisfying condition (A15) . (Cf. the remark made after lemma ^.) 

Assuming that theorems 1 to 3 hold, it is easy to establish 
similar propositions for the derivatives of the function w. 

Theorem A . If t(x) and v(x) are continuously differen- 
tiable, then, using the notation of theorems 1 and 2, the partial 
derivatives q>x, and tyy exist and are continuous in D(a,b), 

whereas qpy exists and is continuous in D(a,b), except for those 
values of y for which K(y) is discontinuous. Furthermore 

|cp x | = |ty|< AT' + BN’ 

|^| < T* + |y| IT’ 

|<p | < A 3 T* + A^BN* 

where T' = max |t'| and N* = max |v 1 | . 

Theorem 5 . If t(x) and v(x) possess continuous 
derivatives of the second order, then, using the notation of 
theorems 1, 2, and 4, the partial derivatives cp xx , i|r X y, 

and iy xx exist and are continuous in D(a,b), whereas the 
derivatives <p X y; cpyy, and tyyy exist and are continuous 
in D(a,b), except for those values of y for which K(y) 
or K'(y) is discontinuous. Furthermore 
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Kx| = feyl ^ AT ” + BN ” 
jcp^l = |%| <A2T" + A 2 |y|N" 

fcpy^J 5 |K*(y)| (AT* + |y|N») + a3t" + A 2 BN" 

where T" = max |t" | and H" =• max |v"). 

Theorem 6 . Under the hypotheses of theorems 3 and A the 
partial derivatives of the first order of q^ and i|r n converge 

(uniformly in D(a,b)) toward those of cp and i|r. Under the 

hypotheses of theorems 3 and 5 the partial derivatives of the 

second order of q^ and i|r n (except perhaps <Pn,yy) converge 

(uniformly in D(a,b)) toward those of q> and ijr. The 

same is true for q> n yy provided that K' n (y) K'(y) uniformly 

for 0 > y > -0. 

The proof of theorems 4 and 5 is based on the fact that 
if w = q) + iijf is a Z-monogenic function, then the functions $ 
and £ defined by 


d s z 


for possible discontinuities 
that if the functions $ 


if) d£Z 

then f is a regular solution of equation (Al) . 

Hence theorems 4- and 5 follow immediately from theorems 1 anfl 2 by 
applying the latter theorems to the Cauchy problems with the ini tial 
data t'(x), V'(x) and t"(x), v"(x). The proof of theorem 6 

proceeds along similar lines and may be left to the reader. 


$ + 


if ~ 


possess continuous derivatives (except 
of $y) satisfying equations (AT), so 

and $ are defined by 


A» 


-A 


$ + if = /(* + 
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It is clear that the functions and derivatives whose existence is 
asserted in theorems 1, 2, k, and 5 exist and are continuous (or piece- 
wise continuous) also along the characteristics hounding D(a,b). In 
fact, the functions r(x) and v(x) and their derivatives have been 
assumed to he continuous in the closed, interval (a,b), so that in 
order to establish the continuity of cp, i{r, i}r x , and so forth in the 

closure of D( a,b) it is necessary only to continue t and v in an 
appropriate manner over the interval (a — e, b + e), e > 0, and 
apply the theorems to the domain D( a — e, b + e). Theorems 1 to 6, 
however, remain valid if t and V and their derivatives are merely 
assumed to be continuous and bounded in the open interval (a,b). 

It is rather obvious what modifications would be necessary if r 
were assumed to be piecewise continuous. The function V may be 
assumed to be bounded and measurable without endangering the validity 
of theorems 1 to as will be seen from the proof. 

4. Formal. Powers . 

A significant- class of S-monogenic functions is obtained by 
setting 

Z^ 0) (z) = 1 

i.z(°)(z) = i 


and,- for n = 1, 2, . . . , 

Z< n )(z) = n J Z z( n “ 1 )(C) d s £ 


i*z( n )(z) = n J* i«z( n 1 ^(^) 


These "formal powers" admit the representation (see reference 12) 
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n n 


Z (n) M =£: (v)x I *- v i' , Y (v) (y) 
V=0 


(A17) 


n n 


i*z( n )(z) = 1 > (V)3p" v i v r^ v )(y) 

V =0 


where 


Y (o) 


(y) = ^°\i) 


= 1 


and, for n = 1, 2, . . . , 


Y (2n + 1 )( y) = (2n + 1} jf* yC^)^) dT] 

Y (2n) (y) = 2n^ 7 K(T 1 )Y( 2n - 1 )(Ti) dr] 

I 

^( 2n+ l)(y)= (2n + 1) J* K^Y^^T)) dtjj 


(A18) 


?^ 2n) (y) = 2n J Y Y^- 1 ^) d n 


Note that the imaginary parts of Z^ and i'Z^ possess continuous 
derivatives even if K(y) possesses discontinuities. 

It follows from the preceding formulas that (for a real x) 
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zW(x) = x 11 


i. Z (n)( x ) = ix °- 


This implies 

Lemma 1. If 


IT , 
T (x) = Y a-^ 

3=o 


v(x) = 21 

3=0 J 


then, for every real number a. 


w(z) 




_b__[z( J+1 )(z) - Z^ +l) (a)] 


(A19) 


is a E— monogenic function satisfying the initial condition 


w(x) = 


Formula (A19) is the analog of D'Alembert's solution of the 
equation of the vibrating string. In fact, if K(y) = — A = Constant, 
equation (Al) becomes 




v(|) d£ + It(x) 






= 0 
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and a simple computation shows that in this case 


Z(*) = g |( x + Ay) n + (x — Aj)^j + ~^{( z + A Y) n — ( 2 : — Ay) n J 
i-z( n ) = A J(x + Ay) n - (x - Ay) 13 ] + i J(x + Ay) n + (x - Ay) 3 ^] 


so that the imaginary paid; of equation (A19) is given hy 


If 


= t(x + Ay) + r(x -J fc ) + 1_ P +Ay v(|) d| 
^ 2 A Jx-Ay 


The analog of Bernoulli’s solution is the one obtained hy super- 
imposing particular solutions of the Chaplygin type. (See section 5 of 
the ANALYSIS . ) 

Consider now two functions, %(y) and K 2 (y) , and denote 
functions (Al8) formed with K = % hy Yj/ n ) and Yj^ 11 ) 
where i = 1, 2. Similarly, let the S(Ki) -monogenic formal powers 
he denoted hy ZjL (n) and .i*Zj^ n ) . 

Lemma 2. Set 


a = a(H) = max |Kf(y)| , 0 'J>y JJ’-H, i = 1, 2 
6 = 5(H) = max |%(y) - ^(y)} , 0 >y 2-fl 


There exist constants M n ’ = M n ’(a) (depending only on a 
and n) such that 
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Yl (n) (y) - Y 2 (n) (y) | < 5V Jy| n 
Yi (n) (y) - Y 2 (n) (y)| < | y] n 


for 0 > y > -H. 

The proof" follows easily hy induction on n, using the definition 
of the function Y. 

Lemma 3 « There exist constants = M n (a) (depending only 
on a and n) such that 


Iz-j^Cz) - Zg( n )(z)| < SMn. |zj ^ 


i-Z^^z) - i.Z 2 ^ n ^(z) 


< 5M n |z| n 


for 0 > y > — B. Here a and 5 have the same meaning as in 
lemma 2. 

The proof follows immediately from lemma 2 and from equations (A17) • 


5. Main Lemma 

The lemma to be established in this section accomplishes a two- 
fold purpose. In the first place, it shows that theorem 3 holds 
whenever theorems 1 and 2 hold for the functions K n and K. In the 
second place, it shows that, if theorems 1 and 2 hold for a sequence of 
functions K n possessing a uniform limit K, then theorems 1 and 2 
hold for K. Since every piecewise continuous function is a uniform 
limit of piecewise constant functions (step functions), it follows 
that theorems 1 to 3 are valid if theorems 1 and 2 hold for the case 
when K is a step function. 

Lemma 4 . Hypotheses: (a) For n = 1, 2, ". . . the func— 

tion K n ( y) satisfies conditions (l) to ( 3 ). (h) For each K^y) 
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theorems 1 and 2 hold, as well as similar theorems concerning the 
Cauchy problem with initial data on the line y = c < 0. 

(c) For 0 > y > -H the sequence K^y) converges uniformly to a 
function K(y) satisfying conditions (l) to (3). 

Conclusion: Let t(x) and v(x) be' continuous functions 
defined for a < x < b. It is assumed that |y| < H in D(a,b) . 
Let w n (z) be a z(K n ) —monogenic function satisfying the condition 


v n (x) 



v(i) d! + iT(s:),' a< x< b 


(A20) 


Then, (a) in D(a,b) the sequence w n converges uniformly to a 
(continuous) function w(z); (b) w(z) is a 2(K)-monogenic; 

(c) w(z) satisfies the initial condition (A15); (d) w(z) 
satisfies the inequalities (Al6); and (e) w(z) is the only 
function in D(a,b) satisfying assertions (b) and (c). 

Bemark .— The functions w n are defined, not in D(a,b), but in 
domain D n (a,b) defined in the following proof. In view of 
equation (A22) , one may speak of the convergence of the sequence w n 
in D( a,b) , since each point of D(a,b) belongs to all but a finite 
number of the regions D n (a,b). 


Proof At first it will be assumed that assertion (a) holds. 
Set v n = + iilr n and w = cp + ity. For every closed curve r 

in D(a,b) the relations 


i 

i 


9 dx — Et dy = 0 


t dx + cp dy = 0 


(A21) 


must hold in order that w(z) be z(K)-monogenic. Let D n (a,b) denote 
the set defined with respect to K n (y) in the same way as D(a,b) was 
defined with respect to K(y). Hypothesis (c) implies that 


lim D n (a,b) = D(a,b) 
n— > 00 


(A22) 
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Hence, for sufficiently large values of n the curve T is 
in D n (a,b). Then 



cbt - Vn = 0 

► 

dx + cp n dy = 0 

J 


(A23) 


so that equations (A2l) follow from hypothesis (c) and assertion (a). 
Assertion (c) follows immediately from equation (A20) . Assertion (d) 
is true, for the estimates, expressions (Al6), hold for each 
function v n according to hypothesis ( a) . In order to prove the 

uniform convergence of the sequence w n , choose a positive e and 
determine two polynomials 


t(x) = j! a*x>3 
3=0 


M 

n(x) = 21 

3=0 J 


such that 


jt(x) - t(x)| < 


Jn(x) - v(x)|< 6 l 


(A2k) 


a < x < b J 

(This is possible by virtue of Weierstrass 's approximation theorem.) 
Let W n (z) be the 2 (K n ) -mono genic function satisfying the condition 
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^n( s ) 



(t) d.| + it(x). 


a ^ x ^ b 


4l 

(A25) 


Since theorem 2 is assumed to hold for K = K n , it follows from 
conditions (A20) , (A24), and (A25) that for all points in D n (a,b) 


Jw n (z) -¥ n (z)j < Le, n = ±, 2, 


(A2 6) 


where L is a constant independent of z, n, and e . Denote by 

zj> m ) and i-Z^ 131 ) the 2 (K n ) -monogenic formal powers (formed with 
respect to the function K n ) . Since W n is uniquely determined hy 
condition (A26) , it follows from lemma 1 that 


W (z) = f- z-zjtkz) + i li-C.^t^’tz) - i.^ J+1) (a)] 

1=5 d i=oi + i L 


Hence lemma 3 and the uniform convergence of the sequence K n imply 
the existence of a number Me such that for all values of z in D(a,b) 


*W Z ) “ ¥ m( z )| < 6 > m >M g , P « 1, 2, . . . (A2T) 


By virtue of expressions (A26) and (A27) 

l^WpC 2 ) - ^ m (z)| < (2L + l)e, m >M e , p = 1, 2, . . . 

for all points of the intersection D(a,b) flDm+ptM) OD m (a,b). 

Since L is fixed and e arbitrary and since D n (a,b) — > D(a,b) , 
assertion (a) follows. 

It remains to prove assertion (e). It must be shown that, 
if g(z) is 2 (K) -monogenic and g(x) = 0, a < x < b, g vanishes 
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identically in D(a,b). It is no loss of generality to assume 

that g(z) is continuous in the closure of D(a,b) (since it would be 

sufficient to consider g in D(a + e, b — e), e > 0) . Then 


lim max jg(x + ic)| = 0, (x + ic)eD(a,b) (A28) 

c — >0 


Furthermore, no loss of generality is involved in assuming that the 
real and imaginary parts of g(z) are continuously differentiable (for 
otherwise g(z) could be replaced by its e(K)— integral) . The argument 
leading to the proof of assertions (a) to (d) maybe repeated to 
establish the existence of a E(k)— monogenic function G c (z) satisfying 
the condition 


G c (x + ic) = g(x + ic), (x + ic)eD(a,b) (A29) 

and such that for x + iy = zeD(a,b) and y< c 

G c (z) = lim G Cjn (z) (A30) 

H — ^ 00 

Gc,n( z ) being a z(K n ) -monogenic function with 

^n^ + ic ) “ s( x + ic), (* + ic)eD( a ,b) (A3l) 

Let z Q = xo + iy Q be any point of D(a,b) with y Q < c. It follows 

from hypothesis (b) and theorem 2 that there exists a number L such 
that 

| G c,n( z o)| ^ L feC* + ic)| 
so that by equation (A30) 


KCOI ^ L |s( x + ic)| 
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and by equation (A28) 


G o( z o) - 0 


To show that g = 0 it will suffice to show that G c (z) = g(z) for all 
values of c. Set, c being fixed, 

f(z) = g(z) - G c (z) 

and 


F(z) = $ + if 


92 /»£ 

Jd+ic- M ; (J d+ic 


where d is such that d+ic lies on the characteristic passing 
through (a,0). The function F(z) is z(h) -monogenic and 


F(x + ic) =0, (x + ic)eD(a,b) (A32) 

The function J possesses continuous derivatives fxx. and $ X y and 1 
an at least piecewise continuous derivative fyy It satisfies 
equation (Al) and, in view of equation (A32) , the initial condition 


$(x,c) = f y (x,c) = 0 , (x,c)€D(a,b) (A33) 


But for y < c, equation (Al) is a hyperbolic type and the classical 
theory yields the result j = 0 in D(a,b), for y< c. (The fact 
that K and K' are permitted to have discontinuities does not 
inpair the validity of this argument.) Thus <5=0 and, 
since f = g = G c . This completes the proof. 
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6 . The Step— Function Case 

In this section K(y) is assumed to he a step function (cf. 
figs. 5 and 6; the latter shows the characteristics and the 
domain D(a,b)). 

Set 


K(y) = Jvl-j . 


(A3b) 


where 


0 = Jo > 7 ± > 7 2 
0 <A 1 </^ . . . 


J 


(A35) 


Theorems 1 and 2 will he proved under this assumption. In view of the 
results of the preceding section, this insures the validity of 
theorems 1 to 3 in the general case. 

Lemma_5. If K is given hy equation (A3 1 *-) , a complex-valued 
continuous function w = cp + itjr is Z(K)— monogenic if and only if 
it is continuous and admits the representation 


w = A v fv( x + V) “ 8vC x - A v y )] + 


I [f v (x + A v y) + g^x - A v y)] , y y+1 > 7 - 7 V (A 36 ) 

(Formulas connecting f v with f v+1 are given in the following 
discussion.) 

Proof.— It is easy to see that the condition is always sufficient 
and that it is necessary when ijr possesses continuous derivatives of 
the first order and piecewise continuous derivatives of the second 
order. Hence, if w is S-monogenic, the function 
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f(z) = 


v(£*) d r t 


d 2 z 


admits a representation of the, form of equation (A 36 ) . Since v = f 
the assertion follows. ^ 

Lemma 6 . If K is given hy equation (A3b) , theorem 1 holds. 
Proof.— Set 


v = q) + i'*l r 


cp = A v |f v (x + Ayy) - gy(x - A v y)| 

♦ = f v ( x + \r) + Sv( x -• Ayy), y v _-L - y - y v 


(A37) 


where 


f i (s) =t T(5) £ v(t,) 

gl u) £*»««> 


dl« 


dl« 


(A38) 
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and for v = 1, 2, . 


fy+lU) 


gy+l^^ 


K i+ 4K 5 -( a --^] + 

2^ “ [■ “ (Vl + M^v] 

K 1_ v^ fv t s+(Av+i+Av)yv ] + 

2 + P + ('Vl ~ A v) J vj 


(A39) . 


It is easy to check that v is E-monogenic in D( a, To) and that 
condition (A15) is satisfied. The uniqueness of this solution 
follows from lemma 5- 

Lemma 7 . If K is given by equation (A3 7 ) , theorem 2 holds. 

Proof .— Because of the linearity of the problem and because of the 
special form of inequalities (Al6) it will suffice to prove these 
inequalities under the as suction that at least one of the functions r 
and v vanishes identically. Assume first that v =0. Using the 
representation of the solution given by equations (A37) to (A39 ) } it 
is seen that 


(| f v + lMsv+l|) - (| f vMsy|)> V = 1, 2, . . . 0 


and that 



g-J < T/2 
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bl 


Hence 


l f V+l| S T / 2 

|«v + l | ST / 2 


and by equations (A37) 


Jq> J - AT 
M ^ T 


since for K given by equation (3*0, A(y) > a v for y> y y-1 

Now suppose that r = 0. In this case the solution admits the 
representation 


where 


v(z) = W v (z) + 21 Wj(z), j ± ^ y ^ y 


(A40) 


W v = $ v + ilv 

Wj = cpj + itj 


* (A4l) 


are Z-monogenic functions determined by the initial conditions 
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^(XjO) = J viz) &Z 

i^x.o) = 0 

a < x H. b 

$ v+l( x ^v) = $ vC^7v) 

SyO^Tv) = 0 

Q,y ^ 2C ^ 1>V ^ V = • • • 

9 v ( x ^v) = 0 

^v( z ^v) = $v( x '7v) 

a^ < x < by, V = 1 3 2j . t . 

Here a v or "by denotes the abscissa of the point of intersection of 

the line y = y v with the left or right characteristic, respectively, 
hounding D(a,b). Note that when K is given by equation (3*0 the 
characteristics of equation (Al) are polygonal paths having, 

for > y> yy , the slopes Ay~^ and -Ay - It is easy to 

verify that 
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a^_ — a = — A^y^_ 

a v — a v_i = -^vCyv— 1 ~ ^v) 3 v = 2, 3j • • • 


For 0 > y iT j l3 


^(2,7) = \ 


ix+Ajy 


v( I) d £ + 


i£ 


x-A±y 


v(e) di 


.%(^y) 


x Px+V^y 
2A L t/x— Aj_y 


v( I) d£ 


Hence 


'^Xjy)! - (* - a)H 


so that, in particular. 


|^( a l^l)| - C a l - a)n 


and 


®l,*( x >3l)[ 


Also 


(^(^y)! < |yi w 


(All-2) 


(a43) 


(a44) 


(a45) 


(A46) 


For y v-1 > y > ^ , v = 2, 3 , . . . 
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§ , (x.y) = ® 

v+l v 3 v 


+ | J Q 


x f Z+A v+ l (7 “ 7v) 


®v,x^ ,y v) d ^ + 


1 P^ v+1 (y-y v ) 

2 J a 


®v 


®v,x( 6 *y v ) d ^ 


5 v+ l(x,y) = — f 

2Ay + 1 Jx 


! p+A v+ i(r-7 v ) 




$ (6 ,y ) cL| 

VjX v 


Hence 


K + x (x ' y) l - |* v (v*v)l + ®“ K.iO^v)! 


so that, in particular. 


VlK+i^ v+ i)| ^ |®v(av^v)| + («v+l “ a v ) “K,x ( 6 ,y v ) l 


and 


|*v+l,x( x ^v+l)| 


Also 


|5v+i( 2 ^)| - |y - 7v| |®v,x(^7v)| 


From inequalities (A4-5) and (A49) it follows that 


K,x»-Jv)| V = 1 , 2, 


(Ah7) 


(A48) 


(a49) 


(A 50 ) 


(A5l) 
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Using this and inequalities (AH) and (Aif8), it is seen that 

Mv y v)l ^ («v - a > 

so that by inequalities (A^3) and (A^7) 

|® v (x,y)| < (x - 7y_i > y > y y , v = 1, 2, . . . (A 52 ) 

Also, by inequalities (A46) , (A50) , and (A5l) 

lf v ( x *y)/ ^ |y - 3V-l| N J 7 v_i ^ 7 £ 7 V , v = 1, 2, . . . (A53) 

so that in particular 


K^V 5 ! ~ y v-l| N ' V = 1, 2 , . . . (A5^) 


Using this last inequality and the previous result on the case v = 0, 
it follows that 





1 


i 


> 


(A55) 


♦3 < |yj - yj_ x n, j » 1 , 2 , . 
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Collecting inequalities (A52 ) , (A53 ) > ami (A55 ) > it is clear that 
|cp(x,y)| j(x - a) + |y|A(y)]]J = B(x,y)H 


|lr(x,y)| |y|U 


This completes the proof. 
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Figure 7-— Subsonic flow past a closed body with localized supersonic regions. 



Figure 8.— Flow through Laval nozzle with localized supersonic regions. 


M < 1 



Figure 9 .— Flow past a curved surface with a local supersonic region. 
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